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Abstract 

A symplecticity-preserving RG analysis is carried out to study a 
resonance structure near an elliptic fixed point of a proto-type sym- 
plectic map in two dimensions. Through analyzing fixed points of a 
reduced RG map, a topology of the resonance structure such as a chain 
of resonant islands can be determined analytically. An application of 
this analysis to the Henon map is also presented. 
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1 Introduction 



Since the renormalization group (RG) method was developed as an asymtotic 
method to study a long-time behaviour of a flow in a dynamical system, pQ 
there have been several attempts to extend the method to apply to discrete 
systems. El Particularly, the application to a symplectic map is of great 
importance since a symplectic map is generally derived as a Poincare map of a 
Hamiltonian flow describing a physical system. However, a naive application 
of the RG method fails to describe a long-time behaviour of a system since the 
naive renormalization process does not secure the symplectic smmmetry to 
an RG map 0j . This difficulty generally arises in a discrete system but does 
not occur in the application to a Hamiltonian flow. Therefore, it is urgent to 
develop the RG method for a discrete system which preseves the symplectic 
symmetry. After some attempts to treat special cases a general proce- 

dure has been proposed recently in order to derive a symplecticity-preserving 
(symplectic) RG map [TU]. The newly proposed procedure is based on the 
observation that a naive RG map is not symplectic but a finite truncation of 
the so-called Loiville expansion of a Hamiltonian flow. Then, a symplectic 
RG map is constructed by means of a symplectic integrator of the underlying 
Hamiltonian flow. 

Employing the newly proposed simplecticity-preserving procedure in the 
RG method, we analyze resonance structures such as a chain of resonant 
islands near an elliptic fixed point in a two-dimensional symplectic map. 
A chain of resonant islands is characterized by alternating positioning of 
hyperbolic and elliptic periodic points in the phase space of a symplectic 
map. In our analysis, these periodic points are shown to be correspond to 
the same number of fixed points of an RG map, which simplifies the analysis 
to prove existence of a chain of resonant islands. As an application of the 
present analysis, resonance structures of the Henon map |7| are reproduced 
by means of symplectic RG maps. A similar analysis to the Henon map 
was tried using the RG method jH]. However, the previous study did not 
succeed in reconstruction of resonance structures of the Henon map in terms 
of a symplectic RG map, which would be the most impressive result of an 
application of the RG method. Here, as an example of general results , we 
carry out a comprehensive study of resonance structures of the Henon map 
within the framework of the RG method. 
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2 Symplectic Map near Elliptic Fixed point 



Let us start with the following symplectic map of action- angle type , (x n , y n ) i— > 
[x n +i, Vn+l)' 

y n+l = y n + f(x n ). (1) 

The origin (x n ,y n ) = (0,0) is assumed to be an elliptic fixed point of the 
map (JTJ). Expanding f(x n ) around the fixed point, we have 

x n+l = x n + y n +l 

y n +i = Vn+Yl fl ™<- ( 2 ) 

m=l 

where constants. Since the origin is supposed to be elliptic, there are 

two eigenvalues at the origin: exp(±zu;), where to is real and given in terms 
of the coefficient a± as 

2 cosu — 2 + ai, (3) 
where > a± > —4. Eliminating y n from Eq.(|2*|l , we get 

X n+ i - 2 COS UJX n + X n -! = a 2 x n + a 3 x n H . (4) 

In the case a m = for m > 2, the map (j3J) is called the Henon map [7]. 
If only the coefficient 03 is retained while all other nonlinear coefficients 
vanish, Eq. (jU ) becomes the double- well type map studied in |3j , [12] . The 
resonance structure near an elliptic fixed point is studied by taking Eq.(JU) 
prototype in the following sections. 

3 Perturbative Rrenormalization Analysis of 
Resonance Structure 

Resonance structures in the phase plane of the prototype map (jlj)are analyzed 
by the perturbative RG method. We focus our attention to the case where 
the frequency to is close to one of resonant frequencies 27r//c, k = 3, 4, 5: 

io = 2-K/k + e p 5 } (5) 
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in which e ia a small parameter representing small resonance detuning and 
p{> 0) will be chosen later so that the small resonamce detuning balances 
with a nonlinear term of the leading order. Then, the map (0J reads 

x n+ i— 2 cos(27r/A;)x n +a; n _i = 2{cosc<j— cos(2%/k)}x n +a2X 2 l +a3X 3 l +- ■ ■ , (6) 

where the first term of the RHS of Eq. (jUJ) represents small resonance detun- 
ing. 

Near the origin of the phase plane, assume x n be expanded as 

then the first order perturbation equation gives 

Lx® = z£i - 2x™ cos(2vr/A;) + x^l x = 0, (8) 

from which we have 

X® = Aexp{i(2n/k)n} + Aexp{-i(2n /k)n}, (9) 

where A is a complex constant and A is the complex conjugate to A. The 
ordinary perturbation analysis provides x^ (m > 2) in terms of polynomi- 
als of A and A . Among various polynomials, secular polynomials, of which 
coefficients depend on some powers of n, play a crucial role for a long-time 
behaviour of the system. In our RG method, a special transformation for 
A called an RG transformation is introduced an d an RG map is 

constructed for the renormalized A so that such secularity is removed. Let 
us estimate the order of magnitude of such secular polynomials. We gener- 
ally encounter the following type of secular polynomial known as a term of 
nonlinear frequency shift. 

{dn\A\ 2 + {C 2 n + C 3 n 2 )\A\ 4 + ■ ■ -)Aexp(ium), (10) 

where C m (m = 1, 2, 3) are constants. Thus, the term of nonlinear frequency 
shift is of the order e 3 or higher. 

Since the frequency of the leading order solution © is the resonant fre- 
quency 27r//c, the following secular term appears, in addition to the nonlinear 
frequency shift term (fTUj) . 

C^nA^, (11) 
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where C4 is a constant and k corresponds to the resonant frequency 2ir /k. In 
this paper, we call the secular term (jllj) a resonant secular term. The order 
of the resonant secular term is e^" 1 ^ . Since A is the coefficient of the first 
order solution A exp{—i(27r / k)n} and 

[Aexp{-i(27r/A:)n}] (fc - 1) = A (k - 1] exp{-i(2n(k - l)/k)n} 

= A (fc - 1) exp{i(27r/A;)n}, 

the (k — 1) -th order nonlinear term [Aexpj — i(2ir / k)n}]^ k ~ 1 ^ causes a secular 
solution. For k > 4, the term of nonlinear frequency shift ()10j) is a dominant 
nonlinear term while the resonant secular term (|TT|) yields small but impor- 
tant contribution . Since the resonance detuning e p S should balance with the 
dominant nonlinear term, p should be chosen as p = 2 in this case and for 
k = 4. In the case k = 3 , the resonant secular term dominates over the the 
term of nonlinear frequency shift and p should be chosen as p = 1 so that 
the resonant secular term balances with the resonance detuning term. 
In the following subsections, we carry out the perturbative RG analysis of 
© for resonant frequencies 2n/k{k = 3,4,5). 

3.1 Resonant Frequency : 27r/3 

For k = 3, the resonant frequency becomes 27r/3 and the small resonance 
detuning is chosen to be of the order e as 

to — oo = e5, 

where u = 2n/3. Substituting the expansion (J7J), where the first order 
solution is given by Eq.®, into the map (0), we have the following 
second order equation. 

Lx^ = — 25 sin uuo{A exp(iu} n) + c.c.} + a2{A 2 exp(2iuon) + c.c. + 2|v4| 2 }, 

(12) 

where c.c. stands for the complex conjugate to the preceding term(s). Since 
exp(±2zco> ) = exp(=Fia;o) for ojq = 2n/3, the nonlinear terms proportional to 
A 2 and A 2 and the linear detuning term in Eq. (jl2|) causes secular terms in 
the second order solution x^ : 

xffl = i(5A — A 2 )n exp(iujQn) + c.c. + a 2 o|^4| 2 , (13) 

2sinco>o 
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where 

a 2 ,o = • 14 

1 — COS OJq 

Up to the second order in e, we have 

x n /e = (A + ie(5A — ai 2 )nj exp(iuj n) + c.c. + ea 2i o|^4| 2 , (15) 

where 

a 2 

a 



2 sin uo 

In order to remove secularity in Eq. (jl5|) and derive a symplecticity-preserving 
RG map, we follow the procedure established in [TU] or j^j. First, we intro- 
duce the renormalization transformation A i— > A r 



A n = A + ie(<L4 - — ^ — -4 2 )n, (16) 
2 sin cj 

from which we obtain a naive RG map A n i— > v4 n +i- 

A n+1 = A n + ie(<L4 n - o^) + 0(e 2 ). (17) 

Then, from Eq.([15j). the original variable x n is rewritten in terms of the 
renormalized amplitude A n as 

x n /e = A n exp(iuj n) + c.c. + ea 2) o\A n \ 2 . (18) 

The naive RG map (j!7)l is not symplectic but is made symplectic by means 
of the symplectic integration method. It should be noted that the map 
(JT7jl is the first order truncation of the following Liouville expansion of a 
Hamiltonian flow 

A(t + e) = exp(eC H )A(t) = M + eC H + £ -C 2 H + ■ • • ) A(t), (19) 

with a Hamiltonian 

(A(t), = i(M(f)^(t) - 1 r^ L - {Aitf + A(t) 3 )). (20) 

O Sill LOq 

where A(t + e) = A n+ %, A(t) = A n and the Liouville operator C H is defined, 
in terms of the Poisson bracket {, }, as 

C H A = {A, H}, 

C 2 H A = C h (ChA)={{A,H},B}. 
6 



Therefore, the Hamiltonian ( |2*U|) describes the phase space structure of the 
renormalized amplitude (A n , A n ) within the order of the present approxima- 
tion, regardless of choise of symplectic integrators to the Hamiltonian H . 
The most important points in the the phase space (A n , A n ) are fixed points 
since they ,via Eq.(fT%]). correspond to periodic points in the original phase 
space (x n , x n+ i), which determine resonance structure such as a chain of res- 
onant islands . The fixed points to the flow of H are solutions of the following 
algebraic equations 

dH(A n ,A n ) _ n dH(A n ,A n ) _ n 

dA n " U ' dA n " U - {21} 

In the present case, Eqs. (}2~T]) are reduced to 

5 - - . rexp(-3zfl) = 0, (22) 
2 sin ujq 

where A n = r exp(i8)(r > 0) and we obtain three hyperbolic fixed points 

2.4. r= 2^ 
3 3 a 2 

for S/a 2 > or 

7r 3n 5n — 2<5sincj , N 

^ = 77,^,^, r = °- (24) 

3 3 3 a 2 

for S/a 2 < 0. These fixed points are also those of the naive RG map (|17[) and 
symplectic RG maps derived later. Thus, we find that the resonance structure 
in the original phase plane is characterized by the three hyperbolic periodic 
points with a period 3 instead of a chain of resonant islands. Higher order 
corrections to the dominant Hamiltonian ( l2*Uj) slightly modify the values 
of fixed points and do not change this result as a whole. Two symplectic 
RG maps are derived from the naive RG map fTTj) by choosing approriate 
canonical variables to the Hamiltonian. For canonical varibles {A' n , A'^) where 
A n = A' n + iA£, A' n , A" n e R, we have 

^n+i = A' n + efs(A' n+lJ A'n), 

A» +1 = < + ^ 3 « +1 ,<), (25) 

where f 3 (A',A") = -(5A" + 2aA'A") and g 3 (A',A") = 5A'-a{{A') 2 -{A") 2 ). 
Introducing the action-angle variables (J n ,0 n ) where A n = \fJnexpi9 n , we 
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have another symplectic RG map. 

Jn+l — Jn + £ F 3 [J n+ i, 8 n ), 

9 n + eG 3 {J n+1 ,9 n ), (26) 

where F 3 (J,0) = -2aJ 3/2 sin(30) and G 3 (J,G) = 8 - ot^fj cos{W). The 
former symplectic RG map (|23j) takes a set of explicit diference equations, 
while the latter map (|26|) exactly preserves a symmetry of 2tc/3 rotation 
which the Hamiltonian (J2(Jj) and the naive RG map (|17|) possess. The RG 
map (|26|) in the angle-action variables was first derived in [S]. 



3.2 Resonant Frequency : 27r/4 

For k — 4, the resonant frequency ujq becomes tt/2 and the small resonance 
detuning is chosen to be of the order e 2 as 

u - uj = e 2 5. (27) 

Substituting the expansion (JJJ) into the map © and taking Eq.(|2~7|) into 
account, we have a similar expression to the second order equation as Eq. (jl2j) . 
in which the frequency detuning term shoud be moved to the third order 
equation. Then, in contrast to the previous case, this second order equation 
does not include secular terms and we have the second order solution as a 
non-secular solution 

= a 2j2 A 2 exp(2iuj n) + c.c. + a 2i0 |A| 2 , (28) 

where 

a 2 2 = —, t s r- (29) 

2(cos(2co>o) — coscuq) 

The third order equation is written as 

Lx^ = (-2<5sinu; + (2a 2 (a 2 ,2 + a 2j0 ) + 3a 3 )\A\ 2 ^Aexp(iujQn) 

+ (2a 2 a 2j2 + a 3 )A 3 exp(3itoon) + c.c. (30) 

Since exp(±3iu Q ) = exp(=|=iu;o) for u = ir/2, all forcing terms in Eq.@5UJ) 
contribute to a secular solution. The third order secular solution is given as 

x {3) = i({5 + a 3il \A\ 2 )A + a 3 ,- 3 A 3 )nexp(iuJon) + c.c, (31) 
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where 

a 3 ,i = -tt^~ — (2a 2 (a 2i 2 + a 2j0 ) + 3a 3 ), 
2 sin co> 

03,-3 = -7T-T- — (2a 2 a 2 ,2 + a 3 ). 
2 sin u; 

The secularity in x^ is removed by the following renormalization transfor- 
mation: 

A n = A + ze 2 ((5 + a 3A \A\ 2 )A + a 3i _ 3 A 3 )n, (32) 
which yields a naive RG map. 

A n+1 = A n + ie 2 ((5 + a 3:1 \A n \ 2 )A n + a 3 ^ 3 Af l ). (33) 

The original variable x n is represented in terms of the renormalized amplitude 

A n as 

x n /e = A n exp(iu n) + ea 2 ^A 2 l exp(2«u;on) + c.c. + ea 2) o|^4 n | 2 - (34) 

The naive RG map (J3*3*j) is derived from the Liouville expansion with a Hamil- 
tonian 

H = i(6\A n \ 2 + ^K| 4 + °^f{At + A*)). (35) 

Fixed points in this Hamiltonian flow, which are also fixed points of the naive 
RG map ()33|) , are determined by the similar algebraic equations as those in 
the previous subsection. The fixed points A n = rexp(z#) satisfy 

sin(4fl) = 0, (36) 



,2 



03,1 + 03,-3 C0S(46») 

which give eight fixed points 



> 0, (37) 



o 5 IX 37T 

r = ^ * = ^T> (38) 

r 2 - - 9 - - — — vr — (39) 

r ~ a| + a 3 ' 4' 4 ' 4 4 ' [6J) 



where the following conditions are assumed 

> 0, 

2a 3 a 2 + a 3 



> 0, 72— > 0. (40) 



Four fixed points (|55|) are found to be elliptic, while the other four fixed points 
(|39|) are hyperbolic. Therefore, if the conditions (|4U|) are satisfied, resonance 
structure in this case consists of a chain of four resonant islands. In [TU] . 
such a chain of resonant islands is shown for the case 02 = 0. Following the 
same procedure as in the previous subsection, we derive two symplectic RG 
maps. 

^n+i = A' n + e 2 f A (A' n+1 , A'n), 

K + i = K + e%(K + i,<), (41) 

where f 4 (A', A") = -(5 A" + (o 3 ,i + a 3 ,_ 3 )(A") 3 + (o 3 ,i - 3a 3 ^)(A'fA") and 
9a{A',A") = 6A'+(a 3A + a 3 ^ 3 )(A') 3 + (a 3 , 1 -3a 3 ^ 3 )(A") 2 A>. This symplectic 
RG map takes a set of explicit difference eqations in the case a 2 = [TU] . 
Another one preserves a symmery of 27r/4 rotation exactly. 

■Jn+1 = Jn + £ 2 F4,(J n+ i, 6 n ), 

dn+1 = 9n + S 2 G 4 { J n +1, On), (42) 

where F 4 (J, 9) = 2a 3 _ 3 J 2 sin(46 l ) and G 4 (J, 0) = 5 + a 3j i J + a 3 _ 3 Jcos(46>). 
3.3 Resonant Frequency : 27r/5 

When the resonant frequency ujq is 27r/5, the small resonance detuning is 
chosen to be of the order e 2 . Then, the similar analysis as in previous sections 
yields the following perturbation solution. While xff is given by Eq. (|2^j) . xffl 
is obtained as 

x n = + a 3t i\A\ 2 )Anexp(iuJon)) + a 3t3 A 3 exp(3iui n) + c.c, (43) 
where 

2a 2 a 2 ,2 + a 3 
3,3 2(cos(3c<j ) — cosc<j ) 

and the secular part of xffi is determined by 

Lx^ = ((2a 2 a 3i 3 + 020-2,2 + 3a 3 a 2 ,2 + a 4 )A 4 ) exp(iuj n) 

+2ia 2 (5 + a 3 i\A\ 2 ) A 2 n exp(2iu>Qn) + c.c. + non-secular tern(s^4) 

where a relation exp(±42Co>o) = exp(q=ia;o) has been used. A secular solution 
of Eq. (jUJ) yields the secular part of xff : 

x^ = — ia^A^n exp(iuj n) + ia^ 2 (5 + a 3> i\A\ 2 ) A 2 n exp(2iu n) , (45) 
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where 

2a2<3<3,3 + 0-20-2 2 + ^ a 3 a 2,2 + «4 

"4 = > 

2 sin a>o 

2 

4 ' 2 cos(2u;o) — coscuo 

All secular terms in Eqs. ()43|) and (|45|) are removed by introducing a renor- 
malization transformation. 

A n = A + ie 2 (5 + a^i\A\ 2 )An - ie 3 a A A 4 n. (46) 

Then, we finally obtain a secular-free perturbation solution of the original 
equation (JUJ) up to 0(e 3 ): 

x n = e(A n exp(iuj n) + c.c.) + e 2 ((a 2i 2A 2 l exp(2iuj n) + c.c.) + a 2 fi\A n \ 2 ) 
+e 3 (a 3i3 A 3 n exp(3iuj n) + c.c), (47) 

with a naive RG map for A n 

A n+1 = A n + ie 2 (5 + a 3 , l \A n \ 2 )A n -ie 3 a i A i n (48) 
= (l + e 2 C H + ---)A n . (49) 

A Hamiltonian generating the naive RG map (|48|) is given as 

H(A, A) = i(5\A\ 2 + ^|i|A| 4 - e^(A 5 + A 5 )). (50) 

Although fixed points in the phase plane (A n , A n ) are obtained using the 
naive RG map ()48)1 or the Hamiltonian (|50|). we derive them through a sym- 
plectic RG map. In terms of the action-angle variables ( J n , 9 n ), the map (}4"Hj) 
is easily made symplectic as 

Jn+l = Jn + E 2 F5(J n+ i, 9 n ), 

9 n +l — 9 n + e 2 G 5 (J n+1 ,9 n ), (51) 

where F 5 (J, 9) = -2ea A J 5 / 2 sin(5fl) and G 5 (J, 9) = 5+a 3A J-ea A J 3 / 2 cos(5#). 
Fixed points of the map ([51)1 are given by 

sin(5#) =0, 5 + o 3 ,i J - ea A J 3/2 cos(5#) = 0, (52) 
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from which we have 10 fixed points 



J 



2mrr 



(53) 



J 



^3,1 a 3,l ^3,1 



03,1 a 3,l ^3,1 



5 ' 
(2m + 1)tt 
5 



(54) 



where m = 0, 1, 2, 3, 4 and it is assumed that 




(55) 



The condition (|B3j) simply states that the resonance detuning 5 should be 
compensated by the nonlinear frequency shift. A set of five fixed points (|53|) 
is found to be elliptic while another set (|53|) is hyperbolic. Therefore, if the 
resonance detuning is set so that the condition is satisfied, the resonance 
structure for 27r/5 resonance is composed of a chain of five resonant islands. 
In the next section, the analysis in the present section is applied to study the 
resonance structure of the Henon map. 

4 Application to Henon Map 

A map (jlj with a m = for m > 2 seems to be rather special but has 
been extensively studied as the Henon map|7jjH]. Let us apply results in 
the preceding sections to the Henon map and compare them with numerical 
calculations. 

In the case of 27r/3 resonance near the origin of the phase plane, our 
analysis indicates that there are only three hyperbolic points of period three 
and a chain of resonant island does not exist . Numerical calculations of the 
Henon map in Fig.l supports these results of our analysis. In the inner region 
of a triangle produced by connecting the three hyperbolic periodic points, the 
present analysis gives a good agreement but is not able to describe a chaotic 
sea in the outer region of the triangle. 

The case of 27r/4 resonance is exceptional in the sense that the former 
condition, which guarantees existence of elliptic periodic points, in (J4(Jj) is not 
satisfied since 03 = for the Henon map. Therefore, our analysis indicates 
that there are only four hyperbolic fixed points near the origin of the phase 
plane. Numerical caluculation confirms existence of the four hyperbolic fixed 
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Fig. 1. (a) phase portrait near the 27r/3 resonance obtained from the RG 
map (EHD for e = 0.1 and 5 = 0.8 

(b) phase portrait near the 27r/3 resonance derived by a numerical calculation 
of the Henon map for the same values of parameters. 



points but there also exist four elliptic points far from the origin of the phase 
space as shown in Fig. 2. Our perturbational analysis can not cover a region 
far from the origin. However, if the third order correction to the RG map is 
taken into account, a formal estimation shows that four elliptic fixed points 
exist in a far region, whose distance from the origin is 0(1/ e). 



Fig. 2. phase portrait near the 27r/4 
resonance calculated from the Henon 
map for e = 0.1 and 5 = 0.8 




In the case of 2tt/5 resonance, we have a beautiful agreement of the analy- 
sis and numerical calculations of the Henon map as shown in Fig.3. While the 
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phase portrait of the RG map has a symmetry of 2tt/5 rotation (Fig. 3(a)), the 
phase structure in the Henon map does not have such symmetry (Fig.3(c)). 
This symmetry breaking is realized owing to the higher harmonic terms in 
the solution (}4*T|) as shown in Fig.3(b). 








f 
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Fig. 3. (a) phase portrait of the RG 
map (J51|) near the 27r/5 resonance for 
e = 0.1 and 5 = 0.8, where A' = 
VJ cos(#) and A" = y/1 sin(0) 

(b) phase portrait near the 2tt/5 reso- 
nance in the (x, y) space derived by the 
RG method 

(c) phase portrait near the 2n/5 reso- 
nance calculated from the Henon map 




-0.2 0.2 0.4 



5 Conclusion 

A symplecticity-preserving RG analysis is carried out to study resonance 
structures near an elliptic fixed point of a proto-type symplectic map in two 
dimensions. Symplectic RG maps are constructed as symplectic integrators 
to the underlying Hamilton flows. Fixed points of the RG map correspond 
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to periodic points of the original map, which determine resonance structures 
such as a chain of resonant islands. It should be emphasized that we only need 
to analyze a distribution of fixed points of the RG map in order to determine 
a topological structure of the resonance region. For the resonant frequency 
u = 2n/k (k = 4,5), the RG map has k hyperbolic fixed points and k 
elliptic fixed points around the origin with the same angle spacing, while the 
case k = 3 is special because the RG map has only three hyperbolic fixed 
points. Thus, we conclude that a chain of k resonant islands constitutes 
the resonance structures for k — 4, 5 but such a chain does not exists near 
the origin for k — 3. Although we did not perform analysis for k > 6 
in detail, we are sure that the same conclusion as in the case k = 4, 5 is 
valid in general. In addition to the exisitence of fixed points, the RG map 
has a significant symmetry property that it is invariant with respect to ujq 
rotation around the origin, although the original proto-type symplectic map 
does not posses such a symmetry. This rotational symmetry is broken due 
to the presence of higher hamonics terms and the symmetric RG map can 
reproduce asymmetric resonance structures in the original phase space. 

The above analytical theory for the proto-type symplectic map are checked 
by numerical calculations of the Henon map. The result of numerical cal- 
culations is shown to give a good agreement with the theoretical result, es- 
pecially for the case k = 5. 
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